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ABSTRACT 

We show that the earliest possible reionisation of the Universe is approximately at 
z ~ 13.5 and the optical depth is t ~ 0.17 in the conventionally accepted A cold dark 
matter (ACDM) model with adiabatic fluctuations of a flat spectrum normalised to 
the cosmic microwave background. This is consistent with the reionisation found by 
WMAP (the apparently earlier reionisation epoch of the WMAP is ascribed to the 
adoption of the instantaneous reionisation approximation), i.e., the WMAP result is 
realised only if reionisation of the universe takes place nearly at the maximal efficiency 
in the ACDM model. 

Key words: intergalactic medium - cosmology:theory - cosmology:observation - 
G alaxy : format ion 



1 INTRODUCTION 



The most remarkable discovery from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) Observation is a very 
early reionisation of the universe at z ~ 20 (Bennett et al. 
2003). The optical depth is r = 0.17 ± 0.04 (Kogut et al. 
2003) indicated by the TE cross power correlation function 
for low £ modes, or r = 0.17 ^ ' * from a global fit including 
both TE and TT power spectra (Spergel et al. 2003). 

It is now generally accepted that early reionisation is 
primarily due to UV light from early OB stars in proto- 
galaxies (Couchman & Rees 1986; Tegmark et al. 1994; Cen 
& Ostriker 1992; Fukugita & Kawasaki 1994; see Loeb & 
Barkana 2001 for reviews of the recent progress and further 
references) . Many calculations have been done with increas- 
ingly improved physical approximations (see Cen 2003; Ra- 
zoumov et al. 2002; Ciardi, Ferrara & White 2003 for the 
most recent calculations). Recent calculations are mostly 
based on iV-body simulations to include many physical ef- 
fects. With the cosmological parameters of ACDM models, 
the reionisation epoch is inferred to be z ~ 8 — 11. Although 
these estimates stand for the 'best estimates' with elabo- 
rated simulations, there are some subtleties that may depend 
on specific assumptions, models of physical processes and 
meshes of numerical simulations. In this paper we study the 
earliest possible reionisation epoch in the ACDM scenario 
under the standard assumptions concerning star formation. 
Namely, we consider the case when all energy from baryons 
that form stars are efficiently injected into the intergalactic 
medium. For this purpose it is appropriate to take a ho- 
mogeneous multiphase universe model, where we reduce the 



detailed physics to a few parameters, which are constrained 
by empirical knowledge, or otherwise left as free parameters. 

The cosmological parameters are now well-constrained 
by virtue of WMAP itself, together with advancement in 
optical observations: Ho — 73 km s _ Mpc - , Q m — 0.25, 
JIa = 0.75 (i.e., the flat universe), Qt = 0.043 (Spergel 
et al. 2003). The fluctuations show the spectrum close to 
flat, n = 0.97 ± 0.03, with the amplitude normalisation 
o~s = 0.8 — 0.9. It seems that there is no much freedom 
to modify these parameters. The formation of structure is 
well described by the Press Schechter formalism (Press & 
Schechter 1974). We assume that a constant fraction of 
baryons form stars once the perturbation is collapsed and 
the Jeans and cooling conditions are satisfied. The max- 
imally allowed baryon fraction that goes to stars is con- 
strained by empirical metallicity of the universe. We as- 
sume the Salpeter initial mass function (IMF), which shows 
a rather high weight for high mass stars, relative to Scalo's 
IMF (Scalo 1986), although this choice is not essential to us 
under our empirical constraint from metallicity. We allow 
the fraction of UV that comes out of OB stars in proto- 
galaxies to vary: the extreme case assumes a 100% escaping 
fraction. 

There are many effects that should be taken into ac- 
count for more realistic reionisation history: inhomogeneity 
of star forming clouds, feedback from stars and expanding 
UV shell are obvious examples. The absorption due to galax- 
ies and Lyman alpha clouds are another. These effects all 
contribute to delay the reionisation epoch, so that the ne- 
glect of these effects are justified for the purpose to estimate 
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the earliest reionisation epoch. We discuss, however, these 
physical effects later. 

We find that quasars are important only in the late 
epoch. So we switch off quasars for our calculation of the 
reionisation epoch. 

We follow the formalism given in Fukugita & Kawasaki 
(1994; hereafter FK). We solve the evolution equation for 
thermal history with UV from collapsed objects taken into 
account. We include the effect of heating to the Jeans mass, 
although it is not very important. A full description of the 
basic formalism is found in FK. We give a brief description 
of our calculation and inputs that differ from FK in the next 
section. 



2 CALCULATION 

The number of collapsed objects within the mass range M 
to M + dM at redshift z, N(z, M), is calculated by the Press 
Schechter formula with the threshold mass density S c = 1.68. 
We take a flat (n = 1) spectrum and the transfer function 
given by Bardeen et al. (1986) including baryons. We take 
the high normalisation as = 0.9, allowed by the WMAP 
result (Spergel et al. 2003), but later show how the result 
depends on as- The smearing is made with a top hat win- 
dow function. We take cosmological parameters from the fit 
of WMAPext+2dGRS (with constant power spectrum) in 
Spergel et al. (2003). 

We assume that a constant fraction / of baryons col- 
lapsed into bound objects form stars once the Jeans condi- 
tion 
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where T e is the electron temperature, T 7 is the temperature 
of the cosmic background radiation and fi is the mean molec- 
ular weight, is satisfied, and the molecular cooling time is 
shorter than the dynamical time (Blumenthal et al. 1984), 
i.e., 
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where Y e = lO^fi^Om 2 /! -1 is the fraction of free elec- 
trons. H2 molecules are fragile, and this condition may have 
to be replaced with the atomic cooling condition T v i r > 10 4 K 
(T v i r : virial temperature) for a more realistic calculation 
(Stecher & Williams 1967; Haiman et al. 1997). The atomic 
cooling condition gives the lower bound on the galactic mass 
as 



M > M a , 



2.1 x lO M M (n m /i 2 )- 1/2 (l + z) 



-3/2 



(3) 



but the results using this condition differ very little from 
those we obtained by setting the molecular cooling con- 
dition. We take the object with max[M mc ,Mj] < M or 
max[M ac , Mj] < M as being collapsed. The cooling time be- 
comes longer than the Hubble time for low redshift 1 + z < 
7(Q. m h 2 ) 1 , where Compton cooling is not efficient, and 
galaxies do not form if 



Z being the metallicity (Blumenthal et al. 1984). The con- 
stant fraction / may be estimated by the balance of the local 
infall rate and the cooling rate, which, however, may depend 
on the details of calculations (e.g., Cen & Ostriker 1992). 
Here we simply constrain the / parameter from the heavy 
element abundance to avoid uncertainties from the use of 
specific models. [Note that our / is the fraction of baryons 
that form stars against collapsed baryons that satisfy cool- 
ing conditions, in contrast to the definition often taken in 
the literature that a similar number is defined by the ratio 
of baryons that satisfy the cooling conditions (which equal 
baryons that form stars) against those in collapsed objects.] 
We assume that stars form according to the Salpeter 
IMF, <f>(M a ). The temperature-mass relation (Bond, Carr & 
Hogan 1986) we adopt is 



T S {M S ) = 6 x 10 4 K min 
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for population II stars (Z < 0.001). The coefficient is re- 
placed with 4.3 x 10 4 K for population I stars (Z > 0.01). 
For 0.001 < Z < 0.01 we consider a mixture of the two popu- 
lations to adjust the mean metalliciy. In principle, somewhat 
higher temperature is expected for so-called population III 
stars, i.e., stars from pristine gas, but the local environment 
of star forming region is quickly enriched once the first stars 
formed: so it is more appropriate to consider population II 
stars here 1 . The fraction of radiation energy produced, e s M s , 
and the time in the main-sequence, t-MS, are , respectively, 
(Bond, Carr & Hogan 1986) 
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where X = 0.76 is the hydrogen mass fraction. Since the 
lifetime of massive main-sequence stars with M 3 > 10M©, 
which are relevant to us, is shorter than the cosmic time for 
z < 50, we can assume that ionising photons are produced 
instantaneously upon the formation of stars. The production 
rate of UV photons with energy e 7 is given by 
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where B{e 1 ,T s ) is the blackbody spectrum, normalised to 



J d£jB(e 



1. For simplicity of calculation we use the 



M > 3.OxlO 14 M (l+2) 3/4 (n 6 h 2 ) 3/2 (n m /i 2 )- 5/4 (Z/O.Ol) 3/2 , (4) © 



1 If we take account of population III stars for Z < 10 — 4 , assum- 
ing that star forming regions are not locally enriched after the 
first stars formed, the reionisation epoch is only slightly delayed 
(^reion = 13, t = 0.16, which are comapred to the numbers given 
in section 3 below). This results from the fact that the tempera- 
ture of population III stars is higher (the coefficient of J^J being 
10), but the number of photons is fewer (Carr, Bond & Arnett 
1984), and these two effects counteract for reionisation (the latter 
effect slightly wins); the efficiency of reionisation is more impor- 
tantly controlled by the energy generation rate, which is fixed by 
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derivative of the Press- Schechter function N, whilst being 
aware that it does not give the formation rate of collapsing 
objects when the destruction rate gives a negative contribu- 
tion to the derivative of N. At early epochs, z > 5, however, 
this never happens for galactic scales or smaller scales (Ki- 
tayama & Suto 1996). Thus, this does not introduce appre- 
ciable inaccuracies in the calculation of early reionisation. 

We follow the chemical evolution of the universe. Stars 
with M a > 4Mq eject heavy elements at a fraction 



= 0.5 - (M 3 /6.3M e )~ for 15 < M S /M Q < 100 

= 0.1 for 8 < M 3 /Mq < 15 

= 0.2 for 4 < M s /Mq < 8, (8) 
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Arnett 1984). The evolution of heavy element 
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We take Z < 0.02 as the limit for population I stars, which 
leads to / < 0.55. 

We then solve the coupled evolution equations for ion- 
isation fraction of hydrogen and helium, d[nmi/nn]/dt, 
d[nn c n/riuc]/dt, d[nnein/riHe]/dt, the photon spectrum 
dn~f(ej)/dt and the electron temperature d[kT e nt,/ /j]/dt. We 
take all relevant atomic processes for heating and cooling. 
The explicit form of differential equations and necessary 
atomic physics are detailed in FK. We follow the thermal 
history from z > 1000 including the recombination epoch. 



3 RESULTS AND DISCUSSION 

The results of calculations are shown in Figure 1 (for hy- 
drogen) and in Figure 2 (for helium). The figure shows 
reionisation of hydrogen and excitation of helium I to II 
at Zrcion = 13.5, taking 90% ionisation as the reionisation 
epoch. We define the Thomson optical depth by 



: J dzcor [ n Hii — 7iHii(no reionisation)] 



(10) 



with or the Thomson cross section. Here, we subtract the 
contribution to r from residual ionisation of the standard 
recombination (the second term in the integrand) to evalu- 
ate only the effect of reionisation. We obtain r = 0.172. We 
consider that these z la i on and r correspond to the earliest 
epoch of reionisation and the maximum optical depth we can 
obtain in the ACDM model with the specified cosmological 
parameters and with flat adiabatic fluctuations normalised 
to cosmic microwave background observations. The optical 
depth is consistent with the WMAP result r = 0.17 ± 0.04, 
but the reionisation epoch is later than z rc i on = 17 the 
WMAP team quoted, which assums instantaneous reioni- 
sation (Kogut et al. 2003). [We note that z rc i on ~ 20 often 
quoted by the WMAP team (Bennett et al. 2003; Kogut et 
al. 2003) assumes a rather unusual reionisation history.] The 
difference in z rc ion comes from the fact that stars formed ear- 
lier than the reionisation epoch contribute to r. Although 



the HI fraction drops quite sharply at z r . 



there is a 



non-negligible contribution from early formation of stars at 
z > Zrcion- If we integrate ffllOB to z=13.5, we obtain r = 0.13, 
i.e. At = 0.04 comes from z > 13.5. 
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Figure 1. Ionisation history for hydrogen. The solid (dashed) 
curve represents number fraction of HI (HII). 
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Figure 2. Ionisation history for helium. The solid, dashed and 
long-dashed curves represent number fractions of Hel, Hell and 
Helll, respectively 



Note that first stars formed at z ~ 35, which is visible 
as a turn on of the HII component. This epoch is consis- 
tent among the authors who calculated reionisation with 
the CDM model (see e.g., Loeb & Barkana 2001). 

In our calculation we have ignored many physical ef- 
fects, which always contribute to delay the reionisation. We 
have assumed as a default that all UV photons escape from 
galaxies. This is not realistic, but this fraction is poorly 
known observationally, ranging from >10% (Steidel et al. 
2001) to < 2% (Heckman et al. 2001). Bianchi et al. (2001) 
argued < 20% from the Lyman a cloud abundance. This pa- 
rameter may also depend on the size, and hence redshift, of 
objects we consider. We take this as a free parameter. If we 
assume an escaping fraction of ionising photons Fuv = 10% 
the epoch of reionisation is delayed to z rc i on = 11 and 
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Figure 3. Dependences of the reionisation epoch (indicated by 
solid curves) and the Thomson optical depth (dotted curves) on 
the parameters of our model, (a) Fjjy (b) erg, (c) H m and (d) C^. 

t — 0.13, as we shall show below. We have assumed that 
stars form when the molecular cooling condition is satis- 
fied. If the atomic cooling condition gives a more important 
epoch for first galaxies (Haiman et al. 1997), the epoch of 
reionisation is delayed, but our calculation shows it is only 
by Az ~ 0.006, or At = -0.005. 

We have assumed the Salpeter IMF which gives a rel- 
atively large number of massive stars. The specific choice 
of this IMF, however, results in little uncertainty, because 
we have set a parameter that controls star formation effi- 
ciency / so that it gives empirically observed metallicity. 
Since dominant metal production is due to massive stars, 
the change caused by a different choice of IMF, say Scalo's 
or top-heavy IMF, is compensated by the change of the / 
parameter. 

We note that our results depend not only on / and Fuv, 
but also on other cosmological parameters, Q m and Qf,, and 
in particular on erg. To show sensitivity to various parame- 
ters we show in Figure 3 the dependences of the reionisation 
epoch and the optical depth upon those parameters. The 
dependence on / is identical to that on Fuv- The depen- 
dence of the reionisation epoch on Q m , for example, arises 
from the fact that it controls the formation of small objects 
in early epochs, but on the other hand the decrease of the 
time interval for a unit redshift as increasing Q m suppresses 
the increase of r. The dependence on fi& appears more non- 
trivial: as baryon density increases, recombination rate be- 
comes faster, delaying reionisation. However, this effect is 
overcome by the increase of the baryon density, and r in- 
creases. One may use this figure to compare our calculation 
with those that use different cosmological parameters. 

The other effects that must be taken into account for 
more realistic calculations are feedback from stars to sup- 
press star formation and the absorption of ionising flux by 
hydrogen clouds and the envelope of galaxies. The former 
leads to a milder slope of the luminosity function compared 
with the dark matter mass function (Dekel & Silk 1986; 
White & Frenk 1991; Nagamine et al. 2001). We studied 
the effect by artificially decreasing the net star formation 
rate for smaller objects so that the collapsed baryonic mass 



function obeys the form of an empirical Schechter function 
at z = 0, assuming M/L ^constant. This suppresses the 
early star formation and the ionisation does not take place 
until Zroion ~ 11 and r ~ 0.12. This is probably an extreme 
case because a sharper faint end slope is expected in the 
luminosity function at a higher z. 

The UV absorption by Lyman alpha clouds and galax- 
ies may delay reionisation. We developed a formalism to 
take account of envelope of galaxies and Lyman a clouds 
(Fukugita & Kawasaki, in preparation). For an early epoch 
such as the one we are concerned with, this effect is negligi- 
ble. The absorption due to galaxies becomes sizable only at 
low redshifts (z<5). Adopting a minihalo model for Lyman 
a clouds (Ikeuchi 1986; Rees 1986), we find that absorption 
by Lyman a cloud is negligible for the ionisation epoch, be- 
cause the onset of the ionising radiation decreases rapidly 
the abundance of Lyman a clouds. 

The ionisation by quasars is important only for z < 
3. Recent Sloan Digital Sky Survey (SDSS) observations 
showed that early bright quasars form exponentially in z 
(Fan et al. 2001). We estimate the early quasar popula- 
tion using the quasar luminosity function (Boyle et al. 2000) 
known at lower redshifts, and put [dn 7 /dz] qso with the aid 
of the empirical spectrum (Vanden Berk et al. 2001) in the 
evolution equation. We find that quasars alone reionise the 
universe only at z = 2.6 although the effect becomes recog- 
nisable at z w 10. In any case the contribution of quasars 
is negligible compared to OB stars in the early reionisation 
epoch. The UV flux from quasars is important to sustain a 
high degree of ionisation at lower redshifts, as envisaged by 
the Gunn-Peterson test, when the escaping fraction of UV 
from galaxies is less than 10%. 

In conclusion, the earliest epoch of reionisation we ob- 
tained is z = 13.5 with the cosmological parameters deter- 
mined by WMAP (a value of 14 is possible within errors - for 
example with fl m = 0.3) and the Thomson optical depth is 
0.17. This is consistent with the reionisation WMAP found. 
The difference in z beweeen our value and WMAP's is as- 
cribed to different approximations: WMAP assumed that 
the entire optical depth comes from z < z la i on , whereas the 
contribution earlier than z — z rc i on is non-negligible in our 
calculation. The realistic UV escaping fraction and the feed- 
back of stars may somewhat lower the optical depth and de- 
lay the recombination epoch. The uncertainty in IMF that 
directly controls star formation is effectively absorbed into 
the / parameter which is constrained by metallicity. We re- 
mark that reionisation takes place sharply once the number 
of ionising photons, which are determined by the thermal 
balance, reaches some critical number and so it is difficult 
to sustain a partial ionisation for an extended period, as was 
assumed in Kogut et al. (2003). 

As the final comment our result is consistent with the 
Salpeter IMF case of Ciardi et al. (2003), who carried out 
a numerical simulation for the purpose similar to ours. As 
a check we calculated the case with S7 m = 0.3, og = 0.9 
and Fuv = 0.2, and we obtained the reionisation epoch of 
z = 12.8, which is compared to ~ 12 (see Fig. 3 of Ciardi et 
al.). This small difference is probably ascribed to feedback 
effects of stars, as discussed above. 
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